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Abstract. The first detection of a ^^S-bearing molecule 
in interstellar space is reported. The J=2— 1 and 3—2 tran- 
sitions of C^^S have been observed toward eight Galactic 
molecular hot cores. From a comparison with other opti- 
cally thin isotopic species of CS, the abundance ratio of 
^"S/^'^S is 115(±17). This is smaller than the solar system 
ratio of 200 and supports the idea that ^^S is, unlike the 
other stable sulfur isotopes, a purely secondary nucleus 
that is produced by s-process nucleosynthesis in massive 
stars. 

Key words: Nuclear reactions, nucleosynthesis, abun- 
dances — ISM: abundances — ISM: molecules — Galaxy: 
abundances 
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, !_( ' 1. Introduction 



Sulfur possesses four stable isotopes, ^^S, '^^S, '^^S, and 
■^^S. In the solar system, the abundance ratios are 
95.02:0.75:4.21:0.021 (Anders & Grevesse 1989). For the 
interstellar medium, Chin et al. (1996) determine abun- 
dance ratios of 24.4(±5.0) for ^^S/^^S, and 6.3(±1) for 
^■^S/^^S. While no variation in the ^'^S/^^S isotope ratio 
is found, the ^^S/^^S ratio may increase with galactocen- 
tric radius. In the local interstellar medium (ISM) '^^S, 
^^S and ^''S abundance ratios are not drastically differ- 
ent from those in the solar system. For '^^S, no interstellar 
data have been obtained so far. 

•^^S, ■^^S, and ^*S are mainly primary products of oxy- 
gen burning (Woosley et al. 1973, Chin et al. 1996); how- 
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ever, especially the '^^S yield retains a substantial sensi- 
tivity to the metalhcity (Woosley & Weaver 1995). ^^S is 
thought to be a purely secondary isotope produced by neu- 
tron captures on the primary sulfur seeds during helium 
and carbon burning (Thielemann & Arnett 1985, Langer 
et al. 1989). This isotope is, hence, formed in a process 
significally different from that of the other sulfur or CNO 
isotopes accessible by means of molecular spectroscopy at 
cm- and mm- wavelengths. 

We have observed eight Galactic hot cores in two rota- 
tional lines of C^^S in order to determine the abundance 
of ■^^S relative to other sulfur isotopes and to relate this 
isotopic abundance pattern to that of the local Galaxy 
4.6 Gyrs ago, i.e. the solar system value. 

2. Observations and results 

The J ^ 2-1 lines of "c32s, C^^S, and C^^S 
(95.016711GHz, Lovas 1984), and the J = 3 - 2 tran- 
sitions of C^^S and C^^S (142.522797 GHz, Lovas 1984; 
for the other line frequencies see Lovas 1992) were ob- 
served in January and June 1996 with the 12 m ra- 
diotelescope of the National Radioastronomy Observa- 
tory (NRAO) on Kitt Peak toward W3(0H), Orion-KL 
and W51(M). The same transitions as well as ^^C^^S 
3 — 2 were observed in March and June 1996 toward 
IRAS 15491-5426, IRAS 15520-5234, IRAS 16172-5028, 
NGC6334A, NGC6334B and also Orion-KL using the 
15 m Swedish ESO SubmiUimetre Telescope SEST on La 
Silla. In all cases, weather conditions were good or ex- 
cellent. SIS receivers were used with sideband rejections 
of 20 dB or better. A main-beam brightness temperature 
(Tmb) scale was established by a chopper wheel method. 

12m NRAO observations. The FWHP beamwidth of the 
12 m antenna was 65" for the 2 — 1 transitions, and 41" 
for the 3 — 2 transitions. From maps of continuum sources, 
we estimate the pointing to be correct within 10". The 



2 



R. Mauersberger et al.: Interstellar S: a probe of s-process nucleosynthesis 



RAS1 61 72 



NGC6334A 




-60 -40 
Velocity (km/s) 



-20 20 
Velocity (km/s) 



-20 20 
Velocity (km/s) 



40 60 80 
Velocity (km/s) 



W3(0H 



RAS1 5520 




-60 -40 -20 
Velocity (km/s) 



-20 20 
Velocity (km/s) 



-20 20 
Velocity (km/s) 



-60 -40 -20 
Velocity (km/s) 



-60 -40 -20 
Velocity (km / s) 



Fig. 1. The observed profiles of the lines of S 2—1 (lower frames), C^°S 2—1 (mid frames), and C^°S 3 — 2 (upper frames). 
The velocity is with respect to the Local Standard of Rest (LSR), the intensity scale is Tmb(K)- Also shown are Gaussian fits 
to the lines. For Orion, a simultaneous fit of the C^^S 2 — 1 line with a series of CH3CCCN lines is shown. 



focus was determined after sunset and after sunrise. As 
backends, we employed a filterbank with 2x128 contigu- 
ous channels. Each half of the filterbank was connected 
to one of two orthogonal polarization channels of the re- 
ceiver used. The channel spacing of 500 kHz corresponds 
to a velocity resolution of ^ 1.6 km s^^ for the J = 2 — 1 
transitions and 1.1 kms~^ for the J = 3 — 2 transitions. 
All spectral lines were measured in position switch mode, 
the reference position being at -1-30' in Right Ascension 
relative to the source position; 30s of integration time on 
the reference position were followed by 30 s on the source. 



15 m SEST observations. The J = 2 — 1 and 3 — 2 
lines of each isotopic species were observed simultaneously 
with two single channel SIS receivers. At A = 3 mm, the 
beam has a width of 51", and at 2 mm, the width is 36". 
The pointing was estimated to be accurate within 5". A 
dual beam switch mode with reference positions offset by 
±11'40" in azimuth was employed. The backend consisted 
of two acousto optical spectrometers with 46 kHz chan- 
nel separation. For the presentation of the C'^^S data, the 
resolution was reduced to 0.5 MHz. 



Results. We subtracted linear baselines from all profiles. 
The resulting spectra of C^^S J = 2 - 1 and C^^S 2 - 1 
and 3 — 2 are shown in Fig.|^. The profiles for the C'^^'S 
3-2 and the '^^C^'^S 2-1 and 3-2 transitions have, 
within the noise, the same lineshape as C'^^S 2 — 1, and are 
therefore not displayed. The spectra toward Orion show 
that the calibration scale of the 12 m telescope and the 
SEST are, within the noise, the same. Confidence in our 
results is strengthened by the fact that both telescopes 
show the same spectral features. In the following we will 
use the 12 m data for Orion, which have a better signal- 
to-noise ratio. Spectral features at the frequency of the 
J = 2 — 1 and 3—2 lines of C'^^S have been detected toward 
all sources observed (toward W3(0H) and IRAS 15491 we 
only searched for the C'^^S 3 — 2 line). We have fitted 
Gaussians to all the lines. The results are displayed in 
Table 0. 

3. The identification of interstellar C^^S 

From Fig.|l| it is evident that close to the 2 — 1 line 
frequency of C'^^S there is molecular emission from an- 
other line toward some of the sources, such as Orion, 
IRAS 15520, NGC6334B and possibly toward W51(M). 
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Table 1. Line Parameters from Gaussian fits 



Source^ 


L>GC 

(kpc) 


J - 


- J 


mK km s 


"lsr 

^ — 1 

kms 


kms 


/(C^*S) 
/(C^eS) 




12c32gC 


W3(OH)'' 
Orion-KL'* 


10.4 


3 


- 2 


49(13) 


-46.8 


4.2 


181(49) 






9.0 


2 
3 


- 1 

- 2 


- 83 
233(16) 


8.5 


4.5 


114(50?) 
104(7) 


46(20?) 


3450(1500) 


IRAS 15491 


6.1 


3 


- 2 


146(20) 


-46.6 


5.1 


108(15) 


57(7) 


3020(370) 


IRAS 15520 


6.2 


2 


- 1 


120(30) 


-41.7 


5.4 


128(32) 


68(17) 


3680(920) 






3 


- 2 


243(20) 






133(11) 


77(6) 


4170(320) 


IRAS 16172 


5.6 


2 


— 1 


200(30) 


-52.0 


6.0 


121(18) 


60(9) 


2980(450) 






3 


- 2 


257(30) 






119(14) 


64(7) 


3170(350) 


NGC 6334A 


7.0 


2 


- 1 


216(30) 


-6.7 


4.7 


108(15) 


52(7) 


3120(420) 






3 


- 2 


228(30) 






154(20) 


83(11) 


4990(660) 


NGC 6334B 


7.0 


2 


- 1 


106(30) 


-4.2 


4.5 


163(46) 


76(22) 


4570(1320) 






3 


- 2 


117(20) 






156(26) 


75(13) 


4510(780) 


W51(M)'^ 


6.5 


2 
3 


" 1 
- 2 


214(23) 
290(40) 


57.2 


9.6 


105(11) 
104(15) 


43(5) 


2420(280) 



a) Coordinates used (B1950.0): W3(0H): 2'" 23"" 16™ 8, 5i950 = 61°38'57", Orion-KL: 5'"32™46':"7, -5°24'24", IRAS 15491: 
15'"49™13=0, -54°26'30", IRAS 15520: 15'"52'"0O!l, -52°34'26", IRAS 16172: 16'"17'"13':3, -50°28'14", NGC 6334A: 17'"17'"32':0, 
-35°44'05", NGC6334B: 17''17'"34=0, -35°42'08", and W51(M): 19''21""26"'4, 14°24'42"; b) vlsr and A«i/2 have been fixed 
in the fit using the values determined from C"^'*S 2 — 1. c) see text; d) measured with the 12 m telescope 



The 3 — 2 line appears to be less contaminated except 
toward Orion and W51(M). It is known that chemically 
younger molecular clouds such as Orion tend to have a 
more complex chemistry than older clouds where mainly 
simple molecules are observed (Helmich et al. 1994). This 
can explain the difFences in line blending. 

Close to the C^^S 2 — 1 transition there are several 
K lines of the J = 23 - 22 transition of CH3C3N (methyl 
cyanoacetylene; Lovas 1984). This symmetric top molecule 
has been previously detected toward TMC-1 (Broten et al. 
1984). Assuming that the population of the correspond- 
ing levels can be described by a Boltzmann law with a 
kinetic temperature of 90 K (derived from the very similar 
CH3CN, Cummins et al. 1983) and using the line strengths 
and level energies from Lovas (1984), we made a Gaussian 
fit fixing the estimated relative intensities of the CH3C3N 
lines for Orion, where line contamination is worst. The re- 
sulting estimate for the intensity of the C'^^S 2 — 1 line is 
given m Tableg. 

The only line in the Lovas (1984) line catalog poten- 
tially contaminating the 3 — 2 line of C'^^S is a transition of 
CH2CHCN (vinyl cyanide; v = 142.527 GHz) displaced by 
4.5 MHz (9.4kms~^) from C^^S. Vinyl cyanide is known 
to be abundant in Orion but probably not in massive star 
forming regions with an "older" chemistry. The frequency 
separation is sufficiently high for most of the soures inves- 
tigated to separate it clearly from the C'^^S emission. 

We have further evidence that our identification of 
C'^^S is indeed correct. As it is expected for the opti- 
cal thin case, the (3 — 2)/(2 — 1) line ratios are simi- 
lar for each of the optically thin species, including C'^^S. 
Also there is no large variation of the line ratios between 
different isotopomers of a given transition for different 



sources. We conclude that our identification of C'^^S is 
justified. The line parameters are, however, best deter- 
mined in sources with an "old" chemical composition (i.e. 
less complex molecules), such as IRAS 16172, NGC 6334A, 
and IRAS 15520. 

4. The isotopic abundances of sulfur 

We concentrate here on the ■^^S/'^^S abundance ratio since 
i2Q32g lines tend to be moderately optically thick in 
Galactic star-forming regions (e.g. Frerking et al. 1980, 
Linke & Goldsmith 1980, Chin et al. 1996); ^^C^^S would 
involve a double isotope ratio and C'^'^S is less easy to 
interpret due to its hyperfine structure. There is only a 
slight Galactic gradient in the '^^S/'^'^S abundance ratio 
from the data of Chin (1995) and Chin et al. (1996). 

Since C'^^S and C'^^S are most probably both op- 
tically thin and since both have very similar molecu- 
lar constants and were measured with the same tele- 
scopes we can assume that the ratio of integrated in- 
tensities (see Table |l|) for a given rotational transition 
equals the respective abundance ratio. From the data in 
Table El we obtain a weighted mean abundance ratio of 
34s/3ffS=115(±17). The standard deviation given is the 
square root of the weighted average variance of the data 
(Bevington & Robinson 1992). Since we cannot exclude, 
at the present time, that part of the scatter in the data is 
caused by a Galactic gradient or by less systematic vari- 
ations from source to source, we believe that this refiects 
the uncertainties better than the standard deviation of 
the weighted mean. The unweighted mean is 128(±25). 
The interstellar value of ^^g/^sg is considerably smaller 
than the solar system ratio of 200 (see Fig.^. There is 
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Fig. 2. The derived ^"S/^^S abundance ratios from the 2 - 1 
and 3 — 2 linos as a function of the galactocentric distance. 
The thick and the dashed lines show the weighted mean of the 
measurements and the standard deviation of the data. Also 
shown is the solar system value of 200. 

no systematic difference between the ratios obtained from 
the 2—1 and 3 — 2 lines. Due to the scatter of our data 

and the narrow range in galactocentric distances it is pre- 
mature to discuss a Galactic abundance gradient. Using 

32g ^ J(13C32S) 12c 

and the ^^C/^^C relation with galactocentric radius from 

Wilson & Rood (1994), we obtain as the weighted mean 
value from our data 32s/36s=3280(±760). Also for these 
isotopes, the interstellar ratio is smaller than the solar 
system value (4520). 

According to the metallicity dependent yields of ^''S 
in massive stars (Woosley & Weaver 1995), this isotope, 
which is produced during oxygen burning, has a charac- 
ter intermediate to primary and secondary isotopes. "^^S, 
on the other side, is produced as a purely secondary iso- 
tope during s-proccss nucleosynthesis in massive stars, 
predominantly in helium and carbon burning (Thielemann 
& Arnett 1985, Langer et al. 1989), with a possible (also 
secondary) contribution from AGB stars. Only a small 
fraction of the produced '^^S is destroyed in later burning 
phases and during the supernova explosion (S.E. Woosley, 
priv. comm.). According to the comprehensive massive 
star models of Woosley & Weaver (1995), the ratio of the 
production factors (i.e. output vs. input) of "^^S over ^^S 
in 15 — 25 M0 stellar models drops roughly from 10 to 1 
when the stellar initial metallicity is increased from Zq/IO 
to Zq. 

A larger solar ^''S/'^^S ratio than the corresponding 
present ISM value is thus in agreement with stellar ''^'S 
yields increasing with time and simultaneously constant 
■'"'S yields during Galactic evolution, i.e. a larger ISM value 
4.6 Gyr ago compared to today. In fact, a similar trend is 
found in two other isotope ratios which are clearly primary 



vs. secondary, namely for ^^C/^^C (0: 89, local ISM: 70) 
and i^O/^O (©: 2700, local ISM: 1900; cf. Henkel et al. 
1995). This analogy suggests a positive Galactic ^"'S/^^S 
gradient, similar to those observed for the mentioned car- 
bon and oxygen isotope ratios. This can be confirmed or 
rejected by observations of C'^^'S and C'^^'S in sources cov- 
ering a larger range of galactocentric radii. 

5. Conclusions 

For the first time we have detected a molecule containing 
■^^S in space. C^^S could be measured in two rotational 
transitions and toward eight molecular hot cores. The in- 
terstellar '^^S/^^S abundance ratio of 115 (±17) is lower 
than in the solar system (200). 

This is in agreement with the idea that ^^S is synthe- 
sized as a secondary isotope by s-process nucleosynthesis 
in massive stars. 

Acknowledgements. R.M. and N.L. were supported by Heisen- 
berg fellowships by Deutsche Forschungsgemeinschaft (DFG 
grants Ma 1450/2-1 and La 587/8-2). It is a pleasure to thank 
S.E. Woosley for illuminating discussions. 

References 

Anders, E., Grevesse, N., 1989, Geochim. et Cosmochim. Acta, 
53, 197 

Bevington, P.R., Robinson, D.K., 1992, Data reduction and 
error analysis for the physical sciences, McGraw Hill 

Broten, N.W., MacLeod, J.M., Avery, L.W. et al., 1984, ApJ 
276, L25 

Chin, Y.-N., 1995, Ph.D. Thesis, Bonn University 
Chin, Y.-N., Henkel, C, Whiteoak, J.B., Langer, N., Church- 
well, E.B., 1996, A&A 305, 960 
Cummins, S.E., Green, S., Thaddeus, P., Linke, R.A., 1983, 
ApJ 266, 331 

Frerking, M.A., Wilson, R.W., Linke, R.A., Wannier, P.G., 

1980, ApJ 240, 65 
Hclmich, F.P., Jansen, D.J., dc Graauw, T., Groesbeck, T.D., 

van Dishoock, E.F., 1994, A&A 283, 626 
Henkel, C, Wilson, T.L., Langer, N., Chin, Y.-N., Mauers- 
berger, R., 1994, in: The structure and content of molec- 
ular clouds - 25 years of molecular radioastronomy, eds.: 
T.L. Wilson, K.J. Johnston, Springer Verlag, p. 72 
Langer N., Arcoragi J.-P., Arnould M., 1989, A&A 210, 187 
Linke, R.A., Goldsmith, P.P., 1980, ApJ 235, 437 
Lovas, F.J., 1984, SLAIM Magnetic Tape version T84, priv. 
comm. 

Lovas, F.J., 1992, J. Phys. Chcm. Rcf. Data 21, 181 
Thielemann F.-K., Arnett W.D., 1985, ApJ 295, 604 
Wilson, T.L., Rood, R.T., 1994, ARA&A 32, 191 
Woosley, S.E., Arnett, W.D., Clayton, D.D., 1973, ApJS 26, 
231 

Woosley S.E., Weaver T.A., 1995, ApJS 101, 181 



This article was processed by the author using Springer- Verlag 
Wr^X A&A style file L-AA version 3. 



